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ABSTRACT 

We have developed a relativistic model for pulsar radio emission and polar- 
ization by taking into account of detailed geometry of emission region, rotation 
and modulation. The sparks activity on the polar cap leads to plasma columns in 
the emission region and modulated emission. By considering relativistic plasma 
bunches streaming out along the rotating dipolar field lines as source of curva- 
ture radiation, deduced the polarization state of the radiation field in terms of 
the Stokes parameters. We have simulated a set of typical pulse profiles, and 
analyzed the role of viewing geometry, rotation and modulation on the pulsar 
polarization profiles. Our simulations explain most of the diverse behaviors of 
polarization generally found in pulsar radio profiles. We show that both the 'an- 
tisymmetric' and 'symmetric' types of circular polarization are possible within 
the frame work of curvature radiation. We also show that the 'kinky' nature 
in the polarization position angle traverses might be due to the rotation and 
modulation effects. The phase lag of polarization position angle inflection point 
relative to the phase of core peak also depends up on the rotationally induced 
asymmetry in the curvature of source trajectory and modulation. 

Subject headings: polarization-pulsars: general-radiation mechanisms: 
non-thermal 



1. INTRODUCTION 



Even though pulsars were discovered more than four decades ago, their radio emission 
process is still not completely understood. The very high degree of linear polarization and 
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systematic polarization positio n angle (PPA) swing of pulsar radiation have been naturally 
invoked in curvature radiation ( )Sturrock!ll97ll ; iRuderman fc Sutherlandlll975l ). In the frame 
work of curvature radiation models the radio emission is believed to be emitted by relativistic 
plasma streaming 'force-freely' along the open field lines of the super-strong magnetic field, 
the geometry of which is assumed to be predominantly dipolar. In the non-rotating approx- 
imation, the velocity of relativistic plasma will be parallel to the tangents of the field lines 
to which they are associated with, and hence emitted radiation beamed in the direction of 
field line tangents. In the rotating vector model (RVM), as the pulsar rotates, observer sight 
line encounters different dipolar field lines, and results in the 'S' shaped PPA swing, which is 



more or less deter mined by the geometry of emission region (IRadhakrishnan fc Cookdll969 
Komesaroff Ill970l ). 



Pulsar radio emission is believed to be coming from mainly open dipolar field lines which 
lie within the polar cap region. The shapes of individual pulses indicate that the entire polar 
cap might not be radiating uniformly. T he sub-pulse modulation can be explained on the idea 
of iso lated sparks on the polar cap (e.g.. IRuderman fc Sutherlandlll975l ; ICheng fc Ruderman 
19801 ). Pulsar average profiles resulted from the summation of several hundreds of indi- 
vidual pulses, have well defined shapes and in general they made up of many components. 
Phenomenologically, pulsar emission is recognized as central 'core' emission arising from the 
region nea r to the magn e tic pole and 'cone' emission arising from concentric rings around the 



Mitra fc Rankin 



200 



21 



pole fe.g..lRankinlll983l . I199CH . Il993l : Mitra fc Deshpandelll999l : iGangadhara fc Guptal 12001 



However, there are some contrary arguments that the emission is 



'patchy' (ILyne fc Manchester! Il988f ). 



Since pulsars are fast spinning objects, the rotation effects such as aberration, retar- 
dation and polar cap currents are believed to be strongly influencing their emission. For 
an inertial observer, in addition to intrinsic velocity along the field line tangents, particles 
will have co- rotation velocity component. Therefore, in the inertial observer frame, the net 
velocity of particles will be offset from the field line tangents to which they are associated 
with, and hence the emission will be aberrated in the direction of pulsar rotation. By tak- 
ing into account of rotation, Blaskiewicz, Cordes and Wasserman (1991), hereafter BCW 
(1991), have proposed a relativistic pulsar polarization model. By assuming a constant 
emission altitude r across the whole pulse, they have predicted that the midpoint of the in- 
tensity profile shifts to the earlier phase by ~ t/tlc with respect to fiducial phase, whereas 
PPA inflection point shifts to later phase by ~ 3r/rLc- The parameter ric — c P/2tt is 
the light cylinder radius, and c is the velocity of light and P is the pulsar rotation period. 
Further, they have shown that the intensity on leading side becomes st ronger than that on 
traili ng side due to rotation , which has strong observational support (ILyne fc Manchester 
19881 ). iHibschman fc Arond (120011 ) further improved the relativistic RVM model by tak- 
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ing into account of induced magnetic field due to polar cap curren ts. iDvka (|2008f ) has 
presented a more simple deriva tion in which he has reproduc ed the IBCWI (ll99ll ) predic- 
tion. Following these deductions iThomas fc Gangadharal (120101 ) have estimated the absolute 
emission altitudes of a few pulsars. The asymmetry in the phase location of conal compo- 
nents with respect to core has bee n interpreted in terms of ro t ation effects such as aberra- 



tion and retardation phase sh i fts (iGangadhara fe Gupta 



Dvks. Rudak fc Harding! 12004 ; lGangadharall2005l ; iKrzeszowski et all 120091 ) . 
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By solving the equation of motion, IThomas fc Gangadharal (120071 ) have also predicted 
that the emissions on leading side dominate over t he trailing side due to rotat i on in duced 
asymmetry in the curvature of trajectory of bunches. iDyks. Wright fc Demorest I (120101 ) have 
analyzed the influence of rotation on shape of pulse profiles of millisecond pulsars and identi- 
fied two opposing effects of corotation: (1) the caustic enhancement of the trailing side emis- 
sion due to squeezing into a narrower component; and (2) t he weakening of the trailing side 
caused by the the smaller curvature of source trajectories. IThomas. Gupta fc Gangadhara 
( 120 lOh have discussed the significance of geometric and rotation effects on the pulsar ra- 
dio profiles and interpreted the as ymmetry in the pulsar rad io profiles, particularly "partial 
cones" (which were first termed by lLyne fc Manchesterlll988l ) in terms of the rotation effects. 
These are notable for their highly asymmetric average intensity profiles and PPA traverses, 
wherein one side of a double component conal profiles is either missing or significantly sup- 
pressed, and the PPA inflection point lies well towards the trailing side. 

Among the several theoretical problems related to exploring the pulsar radio emission 
mechanism, the mostly unexplained observational fact is the high degr ee of circular polariza- 
tion a nd its very diverse behavior. By analyzing average pulsar profiles, iRadhakrishnan fc Rankin 
(I1990l ) have identified two types of circular polarization, namely, 'antisymmetric', where the 
circular polarization changes its sense near the center of the pulse profile, and 'symmetric', 
where the circular polarization will have the same sense across the whole pulse profile. In 
the case of pulsars with antisymmetric type, they found a strong correlation between the 
sense reversal of circular polarization and t he PPA swing, an d speculated that it could be a 
geometric property of emission mechanism. lHan et al.l (119981 ) have noticed that the circular 
polarization is common in pulsars but diverse in nature, and even though it is generally 
strongest in the central or 'core' regions, is by no means confined to central regions. They 
found a strong correlation between the sense of circular polarization and the PPA swing in the 
double-conal pulsars, and no correlation between the sense reversal of circular polarization 
near the center of pulse profiles a nd the PPA s wing in the pulsars with antisymmetric type 
of circular polarization. Further, lYou fc Hani ( 120061 ) have reconfirmed these investigations 
with a larger data. 
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There are two probable origins of circular polarizatio n proposed for the pulsar ra 



diation: 


either intrinsic to the emission mechanism (e.g., 


Michel 1987: Gil & Snakowski 


1990alJbl: 


Radhakrishnan & Rankin 


f990; ( 


3angac 


hara 


199' 


7. 


20101) or generated bv the prop- 


agation effects (e.g., 


Cheng & Ruderman 


1979: 


Melrose ' 


2003 


). Gil. Kiiak & Zvcki (1993) 



have modeled the single-pulse polarization characteristics of pulsar radiation and argued 
that the strong sense reversing circ ular polarization is a natural feature of curvature ra- 
diation. I Cheng fc Ruderman! (119791 ) have suggested that conversion of linear polarization 
to circular polarization is possible due to expected asymmetry between the positive and 
negative ch arged components of magneto act i ve pla sma in the far magnetosphere. On the 
other hand, iKazbegi. Machabeli &: Melikidzd (Il99ll ) have argued that the cyclotron insta- 
bility, rather than propagation effect, is responsible for the circular polarizati o n of p ulsar 
radiation. By considering the rotation of magnetosphere, iLyubarskii &: Petroval (119991 ) have 
shown that the induced wave mode coupling in the polarizat ion-limit i ng reg ion can result 
in circular polarization in linearly polarized normal waves. iMelrosd (120031 ) reviewed the 
properties of intrinsic circular polarization, and the circular polarization due to cyclotron 
instability, and discussed the circular polarization due to propagation effects in an inhomo- 
geneous birefringent plasma. 



Recently, iGangadharal (120101 ) has developed a curvature radiation model by considering 
the detailed geometry of emission region in the non-rotating pulsar approximation. His 
results supports the fact that the circular polarization survives only when there is a sufficient 
gradient in the sub-pulse modulation, and the antisymmetric circular polarization is an 
intrinsic nature of the curvature radiation. He has confirmed the Radhakrishnan fc Rankin 



( 119901 ) correlation between the sense reversal of circular polarization and the PPA swing. 
He has also shown that the sense reversal of circular polarization is by no means confined to 
central core regions. 

Although an extensive pulsar polarimetric studies are available, the pulsar emission 
and polarization are not well understood due to their diverse nature. Despite the rotation 
effects such as aberration and retardation are strongly believed to influence the pulsar radio 
profiles, a complete polarization model including rotation effec ts has not been attempted 



i n the literature. In the work of Thomas and co-authors, and iDyks. Wright &: Demorest 



( 120 101 ) modeled only the total intensity, an d left the po la rization part untouched. O n the 
other hand, relativistic models proposed by iBCWl (jl99ll ). iHibschman fc Arond (120011 ). and 
Dykd (120081 ) deals with PPA only, and not even the linear polarization. Further, in their 
models, emission from the points at which bunch velocity exactly aligns with the sight line is 
only considered, and the emissions from the neighboring points at which velocity lies within 
~ I/7 with respect to sight line are not considered. But such emissions does have influences 
on PPA swing if modulated as we show in our model. 
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For the first time, we have developed a complete polarization model by taking into 
account of a detailed geometry of emission region, rotation effects and modulation. We adopt 
the features of curvature radiation, and incorporated the Gaussian sub-pulse modulation. 
At any instant of time, observer tends to receive the incoherent curvature radiation from a 
modulated beaming region, which constitutes a small flux tube of dipolar field lines. Based 
on simulated profiles, we discuss the combined effect of rotation and modulation on the 
typical pulsar radio profiles. We ascribe the asymmetry in the pulsar radio profiles in terms 
of combined effects of viewing geometry, rotation, and modulation. In § 2 we derive the 
expressions for radiation electric field in frequency domain and the Stokes parameters. In 
§ 3 we present the simulation of typical pulse profiles, and the discussion in § 4 and conclusion 
in 8 5. 



2. POLARIZATION STATE OF THE RADIATION FIELD 

Relativistic plasma streaming 'force-freely' along the super-strong dipolar magnetic field 
lines emit beamed curvature radiation. The curvature radiation model requires an efficient 
plasma bunching to account for the very high brightness temperature of the pulsar radio 
emission, wherein the pl asma bunches of size less than or equal to the radiation wavelength 
can e mit coherently (e.g.. lSturrocklll971tlRuderman &: Sutherlandlll975l ; ICheng &: Ruderman 
19801 ). In our model, we treat the plasma bunch as a single particle of charge q = Ne, where 
e is the electronic charge, and N is the number of particles. In this paper we alternatively 
use source, particle, or plasma bunch as source of radio emission but they all mean the same. 
The emissions from different bunches become incoherent as such emissions do not bear any 
phase relation. 

Consider an inclined and rotating magnetic dipole in the inertial observer's frame (IOF), 
a stationary Cartesian coordinate system-XYZ with neutron star center O as the origin as 
shown in Figure [TJ The angular velocity ft is considered along positive Z-axis, and the 
magnetic axis rh is inclined by an angle a with respect to O. Consider a radiation source S 
constrained to move along the rotating field line f. The velocity v of the source is given by 

v = neb + O x r , (1) 



where b is unit tangent vector to the field line and r is the position vector of the source, 
and their expressions are given in iGangadharal ( 120101 ). The parameter k specifies the speed 
of the source along the field line as a fraction of the speed of light c. The first term on the 
r.h.s. of Equation ([I]) is the velocity in the corotating frame and is in the direction of the 
associated field line tangent, and the second term is corotation velocity. Hence the velocity of 
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the source is offset from the field line tangent to which it is associated with, and is aberrated 
in the direction of pulsar rotation. The parameter k can be deduced from the Equation (JT]) 
by assuming |v| = (3c, 




sin 2 9' sin 2 6 



Qr 



sin 9' cos G 



(2) 



where f3 = a/1 — I/7 2 , 7 is the Lorentz factor of the source, 9' is the angle between r and 
fi, and is the angle be t ween b and the rotation direction e. The expressions for 9' and 
are given in IGangadharal ( 2005 ). 



The acceleration a = dv/dt of source in IOF is given by 

(kc) 2 db kc 2 8k~ _ f , ^ , n 

+ —^b + 2nc{nxb) + ftx(nxr) 



Ibl 89 Ibl 09 



(3) 



where we have used the expression of arc length of the field line ds — \h\d9 — kc dt wherein 
|b| is the magnit ude of the fie l d line tangent and 9 is magnetic colatitude. The expression 
for |b| is given in IGangadhara (|2004h . 



The first term on r.h.s. of Equation ([3]) is the acceleration of bunch due to curvature of 
dipolar magnetic fi eld line, and this is the only term which exists in the absence of rotation 
(IGangadharal 120101 ). The second term is due to a small change in the speed of bunch due to 
motion along field line. The third and last terms are the accelerations due to Coriolis and 
the Centrifugal forces, respectively. 

As the relativistic source accelerates along the rotating field line, it em its curvature 
radia tion whose spectral distribution at the observation point Q is given by (IGangadhara 
2010h 

E(r , u) = 1 s^am r , (4) 



V27T Ro c J -00 £ 2 

where h = {sin^, 0, cos£} is the observer's sight line and ( = a + a with a being the sight 
line impact angle. The parameters (3 — v/c and — a/c are the velocity and acceleration 
of source, respectively. £ = 1 — /3 • n, and R is the distance from the neutron star center to 
observer. 

At any rotation phase 0' = <fi' m of the magnetic axis and for a given emission altitude 
r, the viewing geometry allows observer to receive the beamed emission only from a spe- 
cific region of the pulsar magnetosphere. Observer receives the maximum radiation from 
the emission point at which v and n are exactly aligned, and we define the corresponding 
magnetic colatitude 9 = 9n and azimu t h 4> = <pQ. By considering a slowly rotating (or non- 
rotating) magnetosphere IGangadharal (120041 ) has derived the expressions for 9q and 0o as 
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funct ions of rotation phase <ft' = <fi' m of the magnetic axis (see Eqns. 9 & 11 in iGangadhara 
2004 ). We solve h ■ v — 1 numerically to find 9q and <J)q, as the exact analytical solutions 
become complicated once the effect of rotation is consid ered. Our numerical algorithm starts 
with #0 and (f> Q obtained from the iGangadharal (120041 ) as the initial guess values, and find 
the refined ones by solving n ■ v — 1. 



The magnetic colatitude 9 and azimuth O as functions of rotation phase <p' m are plotted 
in Figure [2] for both non-rotating (dotted curves) and rotating (solid curves) cases by using 
the parameters a = 30°, P = 1 s, r n = r/r LC = 0.05 and a = ±5° In the rotating case both 
the minimum of 9q and inflection point of <pQ shift to the earlier rotation phase compared to 
those in non-rotating case. The absolute phase shifts of both 9q minimum and the inflection 
point of (po f rom the fidu cial phase (cj)' m = 0) are found to be 2.9°, which is about ~ r n , as 
pr edicted by |BC Wl (ll99ll ). For comparison we have superposed the curves (dashed line) due 



to IBCWI (119911 ) model, and find that their approximated analytical solutions are valid only 
over a smaller rotation phase around 2.9°. 

Even though the observer receives maximum radiation from the emission point (9$, 0o) 
at which v and n are exactly aligned, observer receives a considerable radiation from the 
neighboring emission points too due to the finite width of emission beam. It is about ~ 2/7, 
and the boundary of the emission region centered on n is specified by the condition n ■ v = 
cos(l/7). For the computational purpose we discretized the beaming region into 'beaming 
region points' (BRP), and the coordinates 9 and (ft of each point is specified by 9 e and e , 
respectively. Since the coordinates 6 and <j> are orthogonal, their ranges 9 e ^ min < 9 < 9 e ^ max 
and 4> e ,min < < <fie,max can be used to specify the beaming region boundary, where the 
subscripts (e,min) and (e,max) denote the lower and upper boundaries of the emission region. 
By considering 6 — I/7 and 9 + I/7 as initial guess values for 9 e ^ min and 6 e ^ max at <fi = O) 
we solve n-v = cos(l/7) numerically and find the roots 9 e ^ min and 6 e ^ max . Next for any 9 = 9 e 
within the range between 9 e ^ min and 9 £jmax , we consider O — I/7 and O + 1/7 as the initial 
guess values for <fi e ,min and <fi e ,max, and solve again h ■ v = cos(l/7) numerically to find the 



roots 



and 



'e,max 



Note that beaming regions on leading side become broader compared 



to those on corresp onding; trailing ones, but, they are symmetric in the non-rotating case 
(IGangadharal |2010| ). This is because the sourc e trajectories on leading side get squeezed 



whereas on trailing side they get stretched (e.g.. iThomas fc Gangadharal 120070 



Note that Equation (j4j) is the integration of the electric field of radiation emitted by 
the relativistic source along its trajectory. Observer receives the beamed radiation only for 
a small segment of source trajectory, say between the points Pj and Pf along a rotating 
field line (see Figure [1]). As the source moves from Pi to Pf along the rotating field line, 
time changes from £j to tf, colatitude 9 changes from 9{ to 9f, and the rotation phase 4>' of 
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magnetic axis changes from to (f)'*. From the relation ds = \h\d8 = kc dt, we deduce the 
expression for time t : 

t = [ + K , (5) 



KC 

where K is the integration constant. Using the condition t = <p' m /Q for all 9 = 9 e within 
the beaming region, it follows from Equation (EJ that K = {<f)' m /Q) - (j \b\/(K,c)dO) g=g . 
Therefore we have 



and the corresponding 0' = Qt for the source motion along any given field line within the 
beaming region. Hence the argument of the integrand in Equation (jlj) becomes function of 
9 only: 

= i ^ p im»[(°-W<^mm.^ (7) 



where time t has to be replaced by the expression given in Equation ([6j). 
Let 

A = {A,} = M "xK"-/)x/^l , (8 ) 

KC 

where Aj with i = x, y and are the components of A in the X, Y and Z directions, 
respectively (see Fig. CE]). We series expand the components of A in powers of 9 about 9 e , 
and obtain 

A(9) = a t0 + aa(9 - 9 e ) + a l2 (9 - 9 e f + 0[{9 - 9 e f] , (9) 

where with i = x, y and z, and j = 0, 1 and 2 are the series expansion coefficients. They 
are given by 

a l0 = Ai{9 e ) , an = 4(0 e ) , a l2 = ^(9 e ) , (10) 

where A' { and A" are the respective first and second derivatives of A\ with respect to 9 
evaluated at 9 e . Since the expressions of Ojo, an, and are too big, we have not reproduced 
them here. However one can always reproduce them by differentiating A. 

We set the argument of exponential in Equation ([7]), u {t — n ■ r/c} = C, and series 
expand in powers of 9 about 9 e , and obtain 

C{9) = c + Cl (9 - 9 e ) + c 2 {9 - 9 e f + c 3 (9 - 9 e f + 0[(9 - 9 e ) 4 } , (11) 

where Ck with k = 0, 1, 2 and 3 are the series expansion coefficients. They are given by 

c = C{9 e ) , c x = C\9 e ) , c 2 = \c"{9 e ) , c 3 = \c"\0 e ) , (12) 
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where C, C", and C" are the respective first, second, and third derivatives of C with 
respect to evaluated at e . Since the expressions of c , C\, c 2 and c 3 are too big, we have 
not reproduced them here. 

By substituting Equations © and (TTTj) into Equation ([7]), we obtain the components of 
E(u) : 

/+0O 
(a, + and + a^e^+^+^d-d , (13) 
-oo 

where % = x, y and z. Here = - e and £ = ge^^+^V^^oc). 

Next, by substituting the integral solutions So, Si and S 2 given in Appendix-A, we 
obtain 

Ei(u) = Eo(a i0 So + a a Si + a i2 S 2 ) ■ (14) 

The Stokes parameters: I, Q, U and V have been used as tools to specify the polarization 
state of the radiation field. The net radiation that the observer tends to receive at any 
rotation phase <p' m will be incoherent addition of radiation from bunches from all those 
emission points wh ich lie within the beaming region. Following the definitions given by 



Gangadharal (120101 ) . we estimate the resultant Stokes parameters I s , Q s , U s and V^. 



3. SIMULATION OF PULSE PROFILES 

3.1. Emission From Beaming Region 

By assuming uniform source distribution and using the viewing parameters: a = 45°, 
a = 2°, pulsar rotation period P = 1 s, phase (p' m = 0°, normalized emission height r n = 0.01, 
particle's Lorentz factor 7 = 400 and observation frequency v = 600 MHz, we computed the 
Stokes parameters for the radiation field. The contour plots of /, L and V in [9, 0)-plane are 
given in Figure |3j For comparison we present both the cases: non-rotating - panels (a), (b), 
and (c), and rotating - panels (a'), (b'), and (</). The parameters /, L and V are normalized 
with the corresponding maximum value of I within the beaming region. We find the ratio 
of peak intensity in the rotating case to that in non-rotating is about 2. In all panels, the 
contour levels are marked on the respective contours. Also the contours are gray colored in 
such a way that darker the region lower the corresponding parameter value. 

In the non-rotating case, the trajectory of the sources are same as their associated 
dipolar field lines. Hence the contours of /, and L are symmetric, while those of V are 
antisymmetric with respect to (6, 0o) plane. In the rotating case, the contours of /, L, and 
V gets rotated in (0, </>)-plane due to non-uniform aberration within BRP. 
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Using the parameters a = 30°, a = 5°, r n = 0.05, P = 1 s, 7 = 400 and u = 600 MHz 
at discrete rotation phases of the magnetic axis (p' m = —10°, 0°, and 10°, we simulated 
the emissions from the beaming region and presented the contours of V in Figure HI In 
each panel V is normalized with the corresponding peak intensity I pea k within the beaming 
region. The rotation of circular contour pattern is more towards the inner rotation phases 
(4>' m = 0°) compared to that on outer rotation phases (4>' m = ±10°). This is because, the 
perpendicular distance from the pulsar spin axis to the beaming region decreases towards 
the outer rotation phases, and hence the smaller rotation effects. Further, circular contour 
pattern on the trailing side (<p' m = 10°) gets more rotated compared to that on the leading 
side (<fi' m = —10°). This is due to larger asymmetry in the intrinsic curvature of the field 
lines between the smaller and larger 9 parts within the beaming region on the trailing side 
compared to that on the leading side. 

To describe the behaviors of circular polarization and the swing of polarization position 
angle (PPA), we define the symbols: "+/-" for transition of V s from left-handed (LH) to 
right-handed (RH) and "-/+" for transition from RH to LH. Call the counterclockwise swing 
of PPA [dij) s /d(j)' m > 0) as "ccw" and the clockwise swing (dip s /d4>' m < 0) as "cw". 

3.2. Emission With Uniform Source Distribution 

Assuming uniform source distribution and using the viewing parameters: a = 30°, 
a = ±5°, P = 1 s, r n = 0.05, 7 = 400 and v = 600 MHz, we computed the Stokes 
parameters for the curvature radiation. The profiles of I s , L s , V s and ip s are plotted in 
Figure [5J The parameters in each panel are normalized with the corresponding maximum 
value of I s . In both the cases of cr, I s becomes more stronger on leading side (0' m < 0°) of 
the fiducial phase (</>' m = 0°) than on the trailing side. This is because, due to the rotation 
induced curvature, the source trajectories on leading side become more curved than on the 
trailing side. The dip in the intensity near 0' m = 0° is due larger radius of curvature p 
compared to the other regions. The behavior of L s is similar to I s except it's smaller values 
due to incoherent addition of radiation from bunches. We observe that a small quantity 
of circular polarization survives due to the rotation induced asymmetry. The polarization 
position angle is increasing (ccw) in the case of cr = 5°, whereas it is decreasing (cw) in 
the case of cr = —5°. The PPA inflection point, the phase at which \dip s /d(f)' m \ is maximum 
(indicated by an arrow), is found to be shifted to 4>' m = 8.8° for cr = 5° and to (p' m = 8.6° for 
cr = —5°. These shifts are about 3 r n . 



For comparison we have superposed IBCWI (Il99ll ) PPA curves (dotted curves) on our 
simulated PPA curves. Our simulated PPA profile shapes and the shift of PPA inflection 
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point are found to be in good agreement with BCW (1991) model near the central parts 
(4>' m ~ 0) bu t sligh tly deviated at larger rotation phases due to the a pproximation s made 



in the IBCWI (Il99lh model. Also, note that at any rotation phase (p' m , IBCWI (Il99lh model 
considers only the emission from the central point of the beaming region, whereas we consider 
the emissions from the whole of the beaming region. 



3.3. Emission With Modulation 



In general pulsar average radio profiles consist of many components, which could be due 
to emission from plasma columns that are associated with sparks on polar cap. When sight 
line cuts through such emissio ns, it encounters int ensity pattern, which may be treated as 
appr oximately Gaussian s (e.g.. Kramer et al.lll994j ). We consider time independent modula- 
tion ( lGangadharall2010l ) so that any fluctuation in the intensity strength of individual pulses 
will be smoothed out and hence our simulated profiles are expected to resemble the pulsar 
average profiles. Consider a time independent Gaussian modulation in both the polar and 
azimuthal directions: 



f(0, 



9 p f/al] exp[- 



-(0-0 P ) 2 /^] 



(15) 



where (8 P , (f) p ) define the peak location of the Gaussian and /o is the amplitude. The pa- 
rameters oq = we/(2\/ln2) and o$ = w^/ (2\/ln2), where w# and w^ are the corresponding 
full width at half-maximum (FWHM) of the Gaussian in the two directions. The resultant 
Stokes parameters I s , Q s , U s and V s after taking into account of modulation are given in 
Ganeadharal feoioh . 



3.3.1. Simulation of Core Emission 

To explore the effect of rotation on the central core component of pulsar radio profiles, 
we consider a Gaussian modulation having peak at (6 p , <p p ) = (1°, 0°). We have chosen 
the peak of modulation slightly away from the magnetic axis to have modulation in both 
the 9 and directions. For simulation, we used the viewing parameters: a = 10°, a = 2°, 
P = 1 s, r n = 0.02, 7 = 400 and v = 600 MHz. Since the minimum of the coordinate 
6 is ~ 2/3 a = 1.33° in this case, sight line passes through the emission region where 
modulation strength is slightly below its amplitude. To see the combined effect of rotation 
and modulation, we considered = 0.15 and three cases for <jq = 0.15, 0.004 and 0.002, 
and the simulated polarization profiles are given in Figure El In all the three cases of a$, 
intensity profiles are shifted to earlier phase with respect to the fiducial phase <ft' m = 0° while 
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the polarization position angle profile shifts to the later phase due to effect of rotation. Note 
that in the absence of rotation, the minimum of 9q and the antisymmetric point of 0o (i- e -) 
phase at which O = 0) occur at (p' m = 0° and hence the intensity will peak at cj)' m = 0°. But 
in the rotating case intensity peaks shift to earlier phases by about 1.36°, 1.31° and 1.24° in 
the three cases of ag = 0.15, 0.004 and 0.002, respectively. The phase shift of I s peak is found 
to decrease with decreasing ag. The reasons are, due to aberration, the minimum of 9$ and 
the antisymmetric point of 0o an d hence the peak of modulation (phase at which observer 
encounters maximum modulation strength) shift to earlier phase. But the emission due to 
uniform source distribution becomes more stronger on the leading side of the modulation 
peak compared to that on the trailing side. Therefore, the peak of modulated total intensity 
will be further advanced in phase with respect to the peak of modulation. However this 
extra phase shift of modulated intensity peaks with respect to modulation peaks favor the 
formation of a more symmetric component: larger unmodulated emission on smaller radius 
of curvature side will be less enhanced by weaker modulation whereas smaller unmodulated 
emission on larger radius of curvature side will be more enhanced by stronger modulation. 
Note that if the modulation is more steeper then the intensity profile closely follows the 
modulation. Hence the pulse width and the phase shift of peak of the total intensity will 
decrease as we go from ag = 0.15 to ag = 0.002. 

In the case of ag = 0.15, the circular polarization is antisymmetric and the transition 
is from RH (negative) to LH (positive). The leading negative circular is found to be more 
stronger comp ared to that on the trailing positive circular, whereas they are equal in the non- 



rotating case ( lGangadharall2010l ). The asymmetry in the strengths of negative and positive 
circulars can be explained as follows: due to rotation , the pattern of circular polarization 
(see Figure [3]) gets rotated in the (9, 0)-plane and hence an asymmetry is introduced in 
both 9 and <fi directions. Since we used 9 P = 1° and the minimum of 9 that the observer 
encounters is ~ 1.3°, the modulation in the 9 direction always enhance the emission over 
smaller values of 9 (see Figure [3]) compared to that over larger values of 9 within the beaming 
region. Further, since we used <ft p = 0°, modulation in <p direction selectively enhance the 
emissions over smaller values of \<p\ compared to those over larger values of |0| within the 
beaming region. Since a^ = ag and beaming regions are more extended in <fi compared to 
those in 9, the modulation gradient in cf) dominates over that in 9. Further, the magnitude 
of rotation of circular polarization pattern within the beaming region is more for the inner 
rotation phases (phases closer to (p' m = 0°) as compared to that on outer phases (see Figure 
H]). Hence modulation can selectively enhance the leading negative circular over the trailing 
positive circular. Also, because of above said reasons, the phase location of the sign reversal 
of circular is found to be lagging the phase location of the peak of total intensity by a small 
amount. 
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In the case of ag = 0.004, the trailing positive circular becomes more stronger than the 
leading negative circular, which is opposite to ag = 0.15. Even though the effect of rotation 
on the pattern of circular polarization is same as in the case of ag = 0.15, the modulation 
becomes comparatively stronger in the 9 direction as ag <C a^. Hence, on the leading side, 
the net two dimensional modulation selectively enhances emission over the lower left part 
of the beaming region, whereas on the trailing side, it selectively enhance the emission over 
the lower right part (see Figure H]). Hence the positive circular on the trailing side becomes 
more stronger compared to negative circular on the leading side. Due to above said reasons, 
the phase of sign reversal of circular is found to be leading the phase of total intensity peak 
by a small amount. 

In the extreme case of ag = 0.002, modulation becomes much more stronger in the 9 
direction. Therefore the modulation selectively enhances the emission over the lower part of 
the beaming region through out the pulse window. Hence the circular polarization becomes 
almost positive through out the pulse. Again due to asymmetry in the magnitude of rotation 
of the circular polarization pattern with respect to rotation phase (p' m , the survived positive 
circular is found to be more stronger on the trailing side as compared to that on leading 
side. 

In all the three cases, L s almost follows I s except for its magnitude. Further, when 
V s is weaker L s is found to be little stronger and vice versa. In all the three cases, PPA 
swing is 'ccw' and PPA inflection points (indicated by arrows) are found to be shifted to 
later phases by 3.42°, 2.82° and 1.00°, respectively. The phase shift of the position angle 
inflection point is found to decrease with decreasing ag due to the combined effect of rotation 
and modulation. 



BCWI (119911 ) have predicted that due to aberration both the 9 minimum and the anti- 
symmetric point of shift to earlier phase by ~ r n , (1.14° for r n = 0.02), with respect to 
the fiducial phase 4>' m = 0°. They assumed that the centroid of the intensity profile coincides 
with the 9 minimum and the antisymmetric point of 0. Further, by using the parti cle ac- 



celera tion vector, which reflects direction of the electric field vector in time domain, |BCW 



( 119911 ) have also predicted that the shift of PPA inflection point to later phase by ~ 3 r n , 
(3.44° for r n = 0.02). However, in our model, we estimate the radiation field in frequency 
domain. We consider the effect of rotation along with modulation and a detailed geometry 
of emission region which includes finite beaming regions fro m which the observer can receive 



the considerable radiation, which is not considered in the IBCWI ( 119911 ) model. Hence the 



phase shifts of tot al intensity a nd PPA inflection points can be significantly different from 



those predicted bv iBCWl (I1991f ). 



Note that if one considers the retardation (radiation propagation time delay), the emis- 
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sions from the beaming region at any rotation phase <p' m will be ar rived at later time, i.e., by 
delay St — n ■ r / c, and hence phase delayed by S(p' ret = flSt (e.g., iGangadharal 120051 ) . How- 



ever since we have assumed that emissions from the whole beaming region originate from a 
particular altitude r for a given phase <p' m , they will have roughly the same St ~ r/c (which 
at most differ by ~ lCT 8 s between center to boundary within the beaming region). Hence 
emissions will arrive at the same time t r = t e + St, where t e and t r are the emission and 
reception times of the radiation, respectively. Further, since we consider a constant r across 
the whole pulse, the net emission due to whole beaming region at any phase 4>' m will be time 
delayed by the same St ~ r/c. Hence a constant phase delay of S(j)' ret ~ fir jc = r /tlc = r n 
across the whole pulse is introduced. After taking into account of retardation along with 
aberration, the phase shifts of intensity peak and PPA inflection point for example in the 
case 00 = 0.15 of Figure [6] will be 2.50° and 2.28°, respectively. Since retardation just causes 
the shift of entire aberrated profile by S(p' ret to the earlier phase, we have not reproduced 
simulations by combining retardation along with aberration. However the shapes of the 
aberrated profile will be affected if one considers the varying altitude across the pulse. 

Also note that if one considers the modulation which is broader than the one considered 
in the case of = og = 0.15 of Figure El then the aberration phase shift of the intensity 
peak becomes substantially different from r n . Further V s becomes almost symmetric type 
with negative circular through out the profile. This is because as modulation becomes 
broader pulse also becomes broader. Further in Figure [6], we kept constant and varied ag 
for the three cases. On the other hand if one considers the case in which ag is kept constant 
and 00 is varying from more broader modulation to steeper one, then the behavior of total 
intensity and PPA profiles will be similar to Figure [6j But the evolution of V s will be from 
almost symmetric type to antisymmetric type due to above mentioned reasons. Also, if one 
considers the case in which both the o$ and ag vary by the same amount from more broader 
modulation to steeper one, then the behavior of pulse profiles will be similar to the case 
wherein ag is kept constant and o$ is varying from larger to smaller value. Further if one 
considers the negative a then the polarization profiles behave similar to the positive a except 
for the fact that the polarities of V s and swing of PPA profile will be opposite. 

To see the combined effect of rotation and modulation in the case of sight line passing 
through the other side of the emission region which lies towards the magnetic axis, we 
considered the sight line with a = 1°. For simulations we kept the other parameters the 
same as those of a = 2° case of Figure [61 The simulated polarization profiles for a^ = 0.15 
and the three cases ag = 0.15, 0.006 and 0.003 are given in Figure [3 The phase shifts of the 
total intensity peak in the cases ag = 0.15, 0.006 and 0.003 are found to be 1.44°, 1.46° and 
1.53°, respectively. The phase shift of the intensity peak tends to increase with decreasing 
ag, a behavior opposite to the case of a = 2° (see Figure [6]). This is because, since we have 
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chosen 9 P = 1°, the emission point coordinate 9 will be closer to 9 P at the outer rotation 
phases \<fi' m \ > 0°. But the minimum of the emission point coordinate \<f>\ will be closer to 
chosen <p p = 0° for the inner earlier rotation phases. Therefore, the modulation mapped onto 
a broader pulse phase as <jg decreases. Hence, the phase shift of intensity peak increases as 
we go from oq = 0.15 to 0.003. 

In the case of o$ = 0.15 circular polarization is marginally antisymmetric and the 
transition is from negative to positive. The leading negative circular is found to be much 
more stronger than the trailing positive circular compared to o$ = 0.15 of Figure O Due 
to viewing geometry observer tends to receive radiation from the beaming region whose 
emission points in 9 are always less than 9 P = 1° in this case. Hence the net modulation 
selectively enhances emission over the upper left part of the beaming region (see Figure H]) on 
the leading side, whereas it selectively enhances the emission over the upper right part of the 
beaming region on the trailing side. Hence from Figure H] we see that the leading negative 
circular becomes much more stronger than the trailing positive circular. The phase lag of 
the location of sign reversal of circular with respect to intensity peak is found to be more 
compared to that in the case of a = 2° and cr e = 0.15. In the case of a e = 0.006, the negative 
circular becomes even more stronger compared to positive circular, and in the extreme case 
of gq = 0.003, circular becomes symmetric, i.e., only negative circular survives through out 
the pulse. 

In all the three cases, the position angle swing is 'ccw' and the phase shifts of its inflection 
point for a e = 0.15, 0.006 and 0.003 are found to be 4.80°, 4.87° and 4.97°, respectively. The 
phase shift of the position angle inflection point is found to increase with decreasing <jg unlike 
in the case of a = 2° where it decreases with decreasing oq. The observed opposite trend 
in the phase shift of the PPA inflection point with respect to modulation parameter o$ is 
because of opposite trend in the selective enhancement of the emission over a part of the 
beaming region in the two cases of a. 

Note that so far we have considered the general cases of core modulation where the 
modulation peak is located slightly away from the magnetic axis and modulated in both 9 
and <fi coordinates. The more plausible case for core modulation is a Gaussian whose peak 
is located at the magnetic axis, i.e., at 9 P = 0° and depends only in 9. The corresponding 
modulation function that follows from Equation (115"]) is f(9) = fo exp(—9 2 /<jg). If one sim- 
ulates the pulse profiles with this modulation, V s becomes symmetric type: positive circular 
(LH) for positive a and negative circular (RH) for negative a thro ugh out the pulse d ue to 



the selective enhancement. Note that in the non-rotating model ( iGangadharal 120101 ). the 
modulation in only 9 direction will never give the considerable net circular irrespective of its 
gradient, as it always enhance both the positive and negative circulars by the same amount 
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(see panel (c) in Figure [3]). 

3.3.2. Simulation of Conal Emission 

In this section we consider the combined effect of rotation and modulation on the con- 
centric conal emissions. Consider two Gaussians whose peaks are situated at (0 p , <f> p ) = 
(2°, ± 65°) and rest of the parameters the same as in Figure [6j The modulation peak 
locations are chosen such that the sight line passes through the emission regions with mod- 
ulation strengths below its amplitude and the region of maximum modulation encountered 
by the observer lies towards the meridional plane. The simulated polarization profiles are 
given in Figure [8] for the three cases of oq = 0.1, 0.006 and 0.002. In all the three cases 
leading side component becomes more stronger than the trailing side component. This is 
due to the combined effect of enhancement in the intrinsic unmodulated emission and the 
modulation strength on leading side over the trailing side. However the enhancement due 
to unmodulated emission is more prominent. In the IOF, due to aberration, the plasma 
bunch trajectories become more curved on leading side compared to those on trailing side, 
and hence more emission occurs on leading side. The marginal decrease in the strength of 
trailing side component as we go from the case <jq = 0.1 to 0.002 is due to weaker modulation 
that the sight line encountered on trailing side. Further in all the cases of o~q, the trailing 
side component becomes considerably narrower than the leading side component. This is 
because, even though the modulation has roughly the same steepness on both the leading 
and trailing sides, the radius of curvature becomes more steeper in phase on trailing side 
compared to that on the leading side. Hence it results in a broader component on leading 
side. The phase shifts of the mid point of intensity peaks (the cone centers indicated by 
arrows) to the earlier phase in the cases o~q = 0.1, 0.006, and 0.002 are found to be 0.74°, 
0.94° and 1.10°, respectively. 

In the case of o~g = 0.1, negative circular becomes more stronger compared to posi- 
tive circular on leading side whereas it is vice versa on the trailing side. This is because, 
the magnitude of rotation of circular pattern is more for inner rotation phases while the 
net modulated emission is slightly smaller compared to outer phases and hence selective en- 
hancement of outer side circulars within the intensity components. In the case of oq = 0.006, 
positive circular became more stronger compared to negative circular over the leading side 
component and vice versa over the trailing side component. In the case of oq = 0.1, the 
modulation is broader and hence sight line encounters both the parts of modulation which 
are lying towards and away from the meridional plane. As ag decreases, observer tends to 
selectively encounter the parts of modulation which are closer to meridional plane. This 
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selectively enhances the V s which lies towards <p m = 0°. Hence the inner sides of V s become 
more stronger as compared to outer sides over both the leading and trailing sides of I s . In 
the extreme case of ag = 0.002, V s becomes symmetric type, i.e., only positive on leading 
side and negative on the trailing side. 

In all the cases of ag linear polarization L s profile almost follows the total intensity I s , 
and PPA swing is 'ccw'. The distortions or kinks are due to combined effect of rotation and 
modulation on the emissions over the beaming regions. The phase shift of the PPA inflection 
point has become uncertain as the kinks are affecting the central part of the position angle 
curves. However if one considers a case where central core component lies between two conal 
components then the position angle inflection point can be found without difficulty 

As a next case, we select the modulations with peaks at different azimuthal locations 
(4> p = ±40°) on the same conal ring which is considered in Figure El In this case sight line 
passes through the emission regions with modulation strengths below its amplitude and the 
region of maximum modulation encountered by the observer lies away from the meridional 
plane. By keeping other viewing parameters same as in the case of Figure El we simulated 
the polarization profiles for the cases of ag = 0.1, 0.006 and 0.002 and given in Figured The 
small increase in the strength of the trailing side component as we go from ag — 0.1 to 0.002 
is due to an increase in the modulation strength that the observer encounters on the trailing 
side. The phase shift of the cone centers to the earlier phase in the cases ag = 0.1, 0.006 
and 0.002 are found to be 1.07°, 1.06° and 1.08°, respectively. These shifts are found to be 
almost independent of ag unlike in the previous cases where they considerably dependent 
on the ag. In the case of ag = 0.1, V s is antisymmetric on both the leading and trailing 
sides with the outer circular is more stronger than the inner. This is similar to the cases 
ag = 0.1 in Figure El where the modulation is broader. In the case of a e = 0.006, V s is again 
antisymmetric on both leading and trailing sides but outer circular dominates over the inner. 
This is because, sight line encounters the major part of the modulation, which lies away from 
the meridional plane. Hence there is a selective enhancement of the outer circular relative 
to the inner. In the extreme case of ag = 0.002, V s becomes symmetric on both the sides: 
negative on leading side and vice versa on the trailing. These are of opposite behaviors as 
compared to Figure El In all the three cases of ag, L s almost follows the I s with lower values 
similar to the previous cases, and the PPA swing is 'ccw' and shows the kinky behavior. 

4. DISCUSSION 

Pulsar rotation along with modulation and viewing geometry seems to be greatly influ- 
encing the pulsar radio profiles. Due to rotation the trajectories of sources on leading side 
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becomes more curved compared to those on trailing side, and hence leading side unmodulated 
emission always dominate over the trailing one. If one considers an azimuthally symmetric 
cone modulations then the leading side intensity components become more stronger than the 



corresponding trailing o nes due to the rotationa 



plasma trajectories (e.g ., iThomas fc Gangadhara 



y ind uced asymmetry in the curvature of 



20071 ). In support of this there is a strong 



observational result by iLyne fc Manchester! (119881 ). and our simulations clearly confirm it. 
We also find that the leading side intensity components become wider than the corresponding 
trailing ones due to an asymmetry in the gradient of radius of curvature between leading and 
trailing sides. These findings have an observa tional evidence (jAhmadi &: Gangadharal 120021) 
and a similar behavior has been discussed by iDyks. Wright &: Demorest I (120101 ). Note that 
if we consider higher emission altitude in our simulations then the trailing side component 
gets substantiall y weaker or even vanishing. Hence it could serve as an explanation for the 
"partial cones" (IThomas. Gupta &: Gangadharal |2010| ). 



The fact that the phase shift of the intensity components to earlier rotation phases 
and that of the PPA inflection point to later phase is a natural consequence of effect of 
rotation. The phase shifts of the centroid of pulse and that of the P PA inflection points have 
been predicted to be about r n and 3r n , respectively (jBCWlll99ll ). But in their simplistic 
model considered only the emission from the points at which source velocity vector exactly 
aligns with the observer's sight line. However there is a considerable emission from the other 
points of beaming region, which is influenced by rotation and modulation. As a result the 
shifts of intensity component and PPA inflection point will no longer remain as r n and 3r n , 
respectively. We have shown that, due to pulsar corotation, the pattern of emission within 
the beaming region gets rotated in (9, <f>)— plane, and hence an asymmetry is introduced 
in both the 9 and directions. Due to these asymmetries within the beaming region, 
either the antisymmetric or symmetric type circular polarization become possible depending 
upon the viewing geometry and modulation. If the modulation is more steeper and has 
roughly the same gradient in both 9 and 0, then the antisymmetric circular polarization 
is observed. On the other hand the symmetric type circular polarization is more plausible 
when modulation is broader and has roughly the same gradient in both 9 and 0, and also 
when modulation is more steeper in 9 than in 0. But in literature, circular polarization has 
been modeled only in the non-rotating pulsar approximation. Hence only the antisymmetric 
circular polarization was thought to be a natural feature of curvature radiation, and the 
symmetric type circular polarization was speculated t o be a consequence of propagation 
effect (e.g.. fail. Kiiak & Zvckilll993l : faangadharalhoioh . 



Han et all ( 119981 ) and lYou fc Hani (120061 ) have found that the sign reversal of circular 



polarization is not only associated with the central 'core' region but also found over conal 
components as well as at the intersection of conal components. From our simulations of conal 
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components, it is possible to explain all types of circular polarizatio n sense reve r sals a nd its 
association with either increasing or decreasing PPA of Table 3 in lYou &: Hani ( 120061 ). For 
example, consider the Figures [S] and in the case of erg = 0.1. On the leading side, we get a 
case where the circular polarization changes sign from negative to positive with an increasing 
PPA. Again by considering the <Tg = 0.002 cases of Figures [8] and [9], on leading side we can 
get a case where circular polarization changes sign from positive to negative with increasing 
PPA. 

We do confirm the Radhakrishnan & Rankin's (1990) correlation that the sense rever- 
sal of circular polarization from negative to positive is correlated with 'ccw' PPA swing (or 
increasing PPA) and vice versa, and we argue this as a geometric property o f curv ature 
radiation. In pulsars with 'symmetric' type of circular polarization, lHan et al.l (119981 ) have 
not found any correlation between the sense of circular polarization and the PPA swing. Our 
simulations also indicate that the negative circular polarization can be associated with either 
'ccw' or 'cw' PPA swing depending upon the viewing geometry and modulation locations. 
Simil arly the positive circu lar too can b e asso ciated with 'ccw' or 'cw' PPA swing. How- 



ever, lHan et al.l ( 119981 ) and lYou fc Hani (120061 ) have found that many conal-double pulsars 
show a single handed circular polarization over both the components. The negative circular 
correlated with the 'ccw' PPA swing and the positive circular with the 'cw' PPA swing. 

Our simulations of conal components show that if the sight line is missing the modulation 
peaks and the steepness of modulation in the polar (9) direction is much larger as compared 
to that in the azimuthal (</>) direction, then the circular polarization becomes single handed 
on both the leading and trailing sides but have opposite signs (see for example case <jq = 0.002 
of Figures |8] and [9]) . Note that we considered situations where the leading and trailing side 
modulations symmetrically lie on a cone centered on the magnetic axis. On the other hand 
if one considers a situation where the modulations are asymmetrically located on a cone, 
then the correlation between the sense of circular polarization and the PPA swing in the 
case of conal-double pulsars can be explained. For example, by choosing the locations of 
modulations in the case of oq = 0.002 of Figure[9]at <p p = 40° and —65°, one can get negative 
circulars over both the leading and trailing components (see the cases oq = 0.002 of Figures 
El and |H]), and hence an association of negative circular with the increasing PPA can be 
established. 

The 'kinky' type distortion in PPA profile has been fou nd in some normal pulsars and 
more commonly in millisecond pulsars. iMitra et al.l (120001) have attributed this effect to 
multi polar magnetic field while IMitra fc Seiradakisl (120041 ) have speculated that the aberra- 
tion/ret ara^ionj^ejButo emission can cause the dist orted PPA 
traverses. 



Ramachandran fe Kramerl (120031 ) following iHibschman &: Aronsl (120011 ) have pro 
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posed that the discrete jumps in the PPA profiles are due to magnetospheric return currents. 
However, from our simulations it is clear that even with a constant emission altitude across 
the whole pulse profile the 'kinky' behaviors can be produced in PPA traverses. Due to an 
incoherent addition of radiation field emitted from a beaming region, which is affected by 
rotation and modulation, the distortions in the PPA traverses are introduced. The PPA 
traverse under both the core as well as conal components are found to get distorted. But 
there are observational claim s that central core region is more likely to show RVM distortions 
than the conal regions (e.g., lRankinlll983l Il990l ; iRadhakrishnan fc Rankin! Il990l ). 



In this work we considered a constant emission altitude r across the whole pulse to make 
comparison with the earlier results. However varying emission altitude across the pulse can 
be incorporated. We considered only the rotation and time independent modulation effects 
in our simulations, as we are interested in the polarization properties of curvature radiation 
which is of intrinsic origin. We plan to consider the propagation effects, polar cap currents, 
magnetic field sweep back and higher multi polar components of magnetic field on pulsar 
radio emission in our future works. 



5. CONCLUSION 

By developing a relativistic model for pulsar radio emission we have attempted to explain 
the complete polarization state of the curvature radiation. Our model takes into account of a 
very detailed geometry of emission region, rotation and modulation as detailed in section 3.3, 
which have not been much incorporated in the earlier models. Based on our pulse profile 
simulations, we conclude the following: 

1. The phase shift of intensity components to earlier phase and the PPA inflection point 
to later phase, are strongly influenced by the combined effect of rotation and the 
modulation. 

2. The components on the leading side become stronger and broader than those on the 
trailing side because of rotation. 

3. In an unmodulated emission a small quantity of circular polarization survives due 
to rotationally induced asymmetry, but from the point of view of observations it is 
insignificant. 

4. For the very first time we are able to show that the 'symmetric' type circular polariza- 
tion can be obtained within the frame work of curvature radiation. This result is very 
important from the point of view of emission mechanism. 
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5. Both the types of circular polarization: antisymmetric (+/— or — /+) and symmetric 
(+ or — ), can result any where within the pulse window due to the combined effect of 
rotation, viewing geometry and modulation. This might be responsible for the diverse 
nature of circular polarization. 

6. We argue that pulsar rotation combined with modulation can introduce 'kinky' pat- 
terns into the PPA traverses. 



We thank J. L. Han and Pengfei Wang for stimulating discussions, and anonymous referee 
for useful comments. 
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APPENDIX-A: 



Consider the integral 



By changing the variable of integration d = (x/l) + m, and defining the constants / 
and m = —02/(303), we obtain 



'dx 
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where 2; = -wk= ( ci — ) and C/ = T7=e 
For Im(z) = 0, we know 
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where Ai(z) is an entire Airy function of z with no branch cut discontinuities, and 
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where Ai'(z) is the derivative of the Airy function Ai(z). Therefore, we have 

S = U j ■ 

By differentiating equation (IA-20 with respect to c\, we obtain 



(A-4) 



(A-5) 



S 1 = / 1 ?e <(pl * +Pa * a+C3#,) tZi? 



17 



+00 



x ■ — e v 3 / dx 



u 



/3C3 



Jl 



(-2 



- 3/ 9^? 



Differentiation of equation (1A-2|) with respect to C2 gives 
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Fig. 1. — Geometry of radiation emission from a relativistic source S accelerated along a 
rotating field line f (thick curve) in a stationary inertial frame XYZ with neutron star center 
O as the origin. Q is the rotation axis, rh is the magnetic axis, and b is the field line tangent. 
The rotation direction is e and the net velocity is v. The observer's sight line n lies in the 
fiducial plane (XZ-plane). ZPX, ZML, and XL are the great circles centered on O. The 
observation point Q is at a distance R from the emission point P. 
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Fig. 2. — The magnetic colatitude 6 and azimuth O of the emission point as functions 
of rotation phase <p' m in the non-rotating (dotted curves) and the rotating (solid curves) 
cases. The dashed curves are those due to BCW (1991) model. Chosen P = 1 s, and 
r n — r / r LC — 0.05. The dotted and solid arrows in the upper panels represent the location 
of #o minima while in the lower panels the location of inflection point of 0o i n non-rotating 
and rotating cases, respectively. 
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Fig. 3. — The simulation showing the emission from a beaming region due to an uniform 
distribution of sources in the non-rotating (left column panels) and the rotating (right column 
panels) cases at (p' m = 0°. Chosen a = 45°, a = 2°, r n = 0.01, P = 1 s, 7 = 400, and 
v = 600 MHz. 
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Fig. 4. — The simulation showing the circular polarization V from the beaming region in the 
rotating case at three discrete rotation phases. Chosen a = 30°, a = 5°, r n = 0.05, P = 1 s, 
7 = 400 and u = 600 MHz. 



-28 - 





Fig. 5. — Simulated pulse profiles: chosen P = 1 s, 7 = 400, v = 600 MHz and r n = 0. 05. In 
the iIj. panels: solid curves are due to our simulations and the dotted curves are due to lBCW 
(Il99ll ) model. The arrows in the ip s panels mark the polarization angle inflection point. 
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Fig. 6. — Simulated pulse profiles: Chosen a = 10°, a = 2°, P = 1 s, r n = 0.02, 7 = 400, 
v = 600 MHz, /o = 1, # p = 1°, P = 0° and = 0.15. In each panel I s (solid curves), L s 
(dashed curves) and V s (dotted curves) are normalized with the respective peak intensity. 
The arrows in the ip s panels mark the polarization angle inflection point. 




Fig. 7. — Simulated pulse profiles: Chosen a = 1° and the other parameters same as in 
Figure |6j 
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Fig. 8. — Simulated pulse profiles: The parameters chosen are same as in Figure |6] except 
{0 P , <f)p) = (2°, ±65°) and = 0.1. 
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Fig. 9. — Simulated pulse profiles: The parameters chosen are same as in Figure E\ except 
<\> p = ±40°. 



